Microscopic relaxation mechanisms and linear magnetization dynamics 
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Linear magnetization dynamics in the presense of a ther- 
mal bath is analyzed for two general classes of microscopic 
damping mechanisms. The resulting stochastic differential 
equations are always in the form of a damped harmonic os- 
cillator driven by a thermal field. The damping term con- 
tains both the interaction mechanisms and the symmetry of 
the magnetic system. Back transformation from the oscil- 
lator coordinates to the magnetization variables results in a 
macroscopic tensor form of damping that reflects the system 
anisotropy. Scalar Landau-Lifshitz-Gilbert damping term is 
valid only for systems with axial symmetry. Analysis of FMR 
linewith measurements versus frequency, temperature, and 
fllm thickness in NiFe films shows good agreement with a 
combination of slow-relaxing impurity and magnon-electron 
confluence processes. 

76.50.-l-g, 75.70.Ak, 75.30.Hx 



I. INTRODUCTION 

The study of classical magnetization dynamics of fine 
ferromagnetic particles and thin films is of great interest 
in connection with the rapid progress of nanotechnolo- 
gies and high-density magnetic recording. The most im- 
portant research problems in this field include: (a) the 
problem of magnetization reversal from the former equi- 
librium direction to a new one under an applied reversal 
magnetic field, (b) magnetization escape over a barrier 
from one equilibrium direction to another one under ther- 
mal agitation, and (c) the effect of magnetization thermal 
fluctuations on the noise characteristics of magnetoresis- 
tive devices. All these problems depend on relaxation 
properties of the system, i.e., detailed mechanisms of 
magnetization interaction with other degrees of freedom 
(spin waves, elastic vibrations, conduction electrons, im- 
purities, etc.) of material. Previously substantial effort 
has been focused on the study of mostly phenomenolog- 
ical magnetization relaxation (e.g., ||l|, [2|, [Q. [||, 

i)- . . . 

The aim of this paper is to analyze linear magneti- 
zation dynamics in the presence of a thermal bath for 
two general classes of microscopic damping mechanisms. 
This analysis will be restricted to low-level excitations 
where the magnetization remains close to equilibrium. 
Recently, magnetization relaxation by direct coupling 
processes have been analyzed for the case of a general 
anisotropic magnetic system . The resulting stochastic 



differential equation obtained was in the form of a simple 
damped harmonic oscillator driven by a random thermal 
fluctuating field. The coordinate system of this oscillator 
(or, normal mode) reflected the magnetic anisotropy of 
the system. Although specific damping mechanisms were 
not analyzed, the damping parameter was equivalent to 
the relaxation time. Transformation back to magnetiza- 
tion variables yielded macroscopic dynamics where the 
damping term reflected the system symmetry in a tensor 
form @, @, 0. 

Here this analysis is expanded to examine generally 
both direct and indirect mechanisms for coupling to a 
thermal bath. In all cases a stochastic differential equa- 
tion in the form of a damped harmonic oscillator with 
random forces is obtained (and the corresponding ten- 
sor form of damping, as well). In addition, two specific 
relaxation mechanisms are examined in detail: the di- 
rect process of magnon-electron confluence and the indi- 
rect process of 'slow-relaxing' impurity. It is shown that 
these two processes in combination can explain experi- 
mental measurements of ferromagnetic resonance (FMR) 
linewidth in permalloy thin films. Using reasonable es- 
timates for the parameters associated with these mech- 
anisms, very good fitting is obtained to the linewidth 
dependencies on temperature, frequency and film thick- 
ness. 

A phcnomcnological tool to study magnetization re- 
laxation, which has had widespread use, is the Landau- 
Lifshitz equation |Q] or, its modification with Gilbert 
damping These (LLG) equations conserve the ab- 
solute value of the magnetization (|M| = const) in a 
single domain. They are relatively simple and there- 
fore have been utilized for various calculations and mi- 
cromagnetic simulations. In particular, these equations 
have been used to analyze the frequency dependence of 
FMR linewidth measurements in permalloy films (e.g., 
jm , ) . Although, this data can be fit using the LLG 
equations (with a constant, phenomenological damping 
a), this fact can not be considered as an experimental 
verification for validity of the general nature of the LLG 
approach. As we show in detail, the general tensor form 
of the macroscopic equation reduces to the LLG form 
only for axially symmetric systems. 

The paper is organized as follows. In the Sec. II we 
formulate the basic gyromagnetic motion in a rotating 
frame and demonstrate that gyromagnetic magnetization 
dynamics is equivalent to the motion of a harmonic oscil- 
lator. Physics of a damped oscillator by direct coupling 
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to a thermal bath is reviewed and expanded in Sec. III. 
In Sec. IV the general mechanism of indirect coupling is 
analyzed. An overview of magnetization dynamics with 
a brief comparison to the phenomenological LLG form 
is given in Sec. V. Experimental FMR data for NiFe 
films are analyzed in Sec. VI. It is shown that relaxation 
processes due to magnetization interactions with slow- 
relaxing impurities and conduction electrons fit the data 
well. A discussion in given at the end of the paper. 

II. FORMULATION OF BASIC 
GYROMAGNETIC MOTION IN ROTATING 
FRAME 

A. Small magnetization motion 

Let us consider small-amplitude magnetization mo- 
tions of a single-domain ferromagnetic particle in the 
vicinity of equilibrium state M =MsZo. Here Zq is the 
unit vector in the equilibrium direction and Mg is the 
saturation magnetization. The magnetization rotation 
around effective field in this case, in general, is ellip- 
tical and the magnetic energy can be represented as a 
quadratic form (e.g., [§1, ||]): 

Here Xq and yo ^-re the unit orthogonal vectors in the 
plane perpendicular to the equilibrium direction. H^o 
and Hyg are positive fields, which include both micro- 
scopic and shape anisotropics and the external magnetic 
field and V is the particle volume (see Appendix I for an 
example) . 

B. Normal mode approach 

We begin with a common normal mode approach (e.g., 
psf , ]l^), where the small oscillations of the trans- 
verse magnetization coordinates are first transformed to 
"rotating" complex variables by a linearized Holstein- 
Primakoff transformation [|l5| : 



~ - ^hjM,/2V{a* + a) , (2) 
My,~i^h-fM,/2Via* -a). 

The magnetic energy (|l|) now can be rewritten in the 
quadratic form: 

£/h = Aa*a+iB/2){aa + a*a*), (3) 

where A = j{H,„ + Hy,)/2 and B = 7(il,„ - Hy„)/2. 
This classical transformation yields a form, common to 
quantum mechanics, where a and a* correspond to cre- 
ation and annihilation operators so that in (0) the term 



a* a = 1 — Mzo/Mg. Note that for an isotropic system 
where Hx^ = Hy^, the term B vanishes. 

The dynamic precession equations are given by 

da/dt = -iAa - iBa* , da* / dt = iAa* + iBa. (4) 

The mixed terms in (^ can be eliminated by the linear 
canonical transformation: 

a = uc + vc* , a* = uc* + vc, (5) 
l A + ujQ __B I A-luq 

The magnetization components in terms of c and c* can 
be expressed as 

c - {hjMs/2Vf^ {HyjHxS'^ (c* + c), (6) 
My„ z {njM,/2V)'^' [HxjHy^f^ (c* ~ c). 

With this transformation the energy in terms of the nor- 
mal mode coordinates is simply: 

£ /fi — ijJqc*c, (7) 

where 

L.0 = ^A^ - 62 ^ iVKJl^o (8) 

is the ferromagnetic resonance (FMR) frequency. The 
dynamic equations in terms of c and c* are now indepen- 
dent: 

dc/dt = —iujQC, dc* / dt — iluqc* . (9) 

In order to derive a damped motion for this oscillator it 
is necessary to consider the interaction with a thermal 
bath. 



III. DIRECT COUPLING TO A THERMAL BATH 

Linear relaxation can be most simply described in the 
thermal bath approximation which represents coupling to 
a reservoir in thermodynamic equilibrium. The thermal 
bath model has the following restriction: the interaction 
with the thermal bath is assumed to be weak. Thus, dy- 
namics of a single-domain ferromagnetic particle in the 
zeroth (no thermal bath) approximation represents the 
gyromagnetic rotation of magnetization around an effec- 
tive field. This means that the symmetry of the gyro- 
magnetic motion is dominant in the system and imposes 
some restrictions for damping and fluctuations. Mathe- 
matical methods for a dynamic system interacting with a 
thermal bath have been developed in detail in quantum 
optics (see, e.g., fill, !§). Here we apply these 

methods for magnetic particle dynamics. 
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A. Direct linear coupling with a thermal bath 

We review 0] and consider classically a general linear 
interaction of a harmonic oscillator with a set of har- 
monic oscillators that represent the thermal bath. The 
Hamiltonian for small oscillations {a*, a} (see, (||)) inter- 
acting with harmonic oscillators {5^, bk} has the form: 

n/n^£/h + Y,^kblbk (10) 
fe 

+ Y,iGk{abl + c.c.) + Fk{abk + c.c.)]. 

k 

The first term is the energy of the magnetic system (|^), 
the second term gives the oscillations of the thermal bath 
with coordinates (modes) bk and 6^, and the last term 
gives a direct interaction coupling. The last term, if 
transformed back to the magnetization variables, can be 
represented as a scalar product M • Hth, where Hth is an 
effective field from the thermal bath. 

Applying the transformation (|^) , we can write the dy- 
namic equations for c and bk as 

dc/dt^ ~iu;oc-iY,iGkbk + Fkbl), (11) 

k 

dbk/dt = —iujkbk — iGkC — iFkC* . (12) 

Here Gk = GkU + Fj-v and Fk — GkV + FkU. The terms 
Fkb^, in ( prj ) and iFkC* in (12) vanish as fast oscillating 



terms (they describe rotations in the complex plane in the 
opposite direction to c and bk , respectively) and therefore 
can be omitted (see, [Q for details). A solution of (|lj) 
is: 

bk{t) = 6fe(0)e-'-'=* - I f Gkcit')e-''^-^'-*'Ut'. (13) 



Substituting (13) into (pT|), yields the integro-differential 
equation for at): 



dc 
di 



Here 



-iLUoc - V IGfcp /* c(t')e-^"'=(*-*')rft' + fit). 

I, "'0 



f(.t)^-iY,GkbkiO)e- 



■ibJkt 



(14) 



(15) 



describes a noise: bk (0) are random (and in thermal equi- 
librium). A solution of (Q) can be represented as c{t) = 
c{t) eyip{—iuJot), where c(i) is a "slow" variable. Neglect- 
ing the memory for the slow variable (c(t') — *■ c{t)), we 
can rewrite Eq.(|l3) as: 



dc 

dt 



-ILUqC 



i. "'O 



(16) 



For a long time t ^ l/(wo — w^) one can put the upper 
integral limit at oo and use the formula: 



duexp(iQu) = ttS{Q) 



.P.v.(l 



(17) 



where 6(^1) is the Dirac delta function and P. v. is the 
principal value of the integral. Thus, we obtain the dy- 
namic equation of a damped harmonic oscillator: 



where 



dc/dt + rjc = —i{u!o + Auj)c + f{t), 



(18) 



(19) 



Aw = -P.v, 



\Gk? ^ 



The damping parameter 77 (|19|) represent a particular 
form of "Fermi's Golden Rule" (e.g., ||). This form, 
for example, has been used for calculating two-magnon 
scattering [E0[, |2^] and magnon-phonon [E2| processes. 



B. Arbitrary coupling to a thermal bath 

Lax ||l^ , ||l^ has developed the most general approach 
for a dynamic description of a harmonic oscillator inter- 
acting directly with a thermal bath (reservoir). He used 
the following Hamiltonian: 



n 



total 



= huji0z+ {cf'^ + c'^F) + nR, 



(20) 



where and c are creation and annihilation Bose opera- 
tors, F and its Hermitian conjugate F'^ describe thermal 
bath variables, T-Lr is the reservoir Hamiltonian. Tak- 
ing the density matrix equation (i.e., one of the most 
fundamental equations in physics) , Lax derived the mas- 
ter equation in the second approximation on interaction 
and finally showed that the damped harmonic oscillator 
should obey the following classical Langevin equation: 



dc/dt + rjc= —i{ujo + Aw)c -I- f{t). 



(21) 



Here c — (c) is the complex amplitude (classical analog of 
annihilation operator) and f{t) is a noise. The damping 
77 and frequency shift Aw are defined by 



77 — iAlu 



due-^"«"([F(0),i?t(u)])«, (22) 



where {■■■)r denotes averaging over reservoir. As before 
the formula (22) for elementary processes can be reduced 
to the "Fermi's Golden Rule" (second order perturbation 
theory, e.g., . ||l|l). 
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1. Magnon- electron scattering 

As an example of magnetization relaxation due to di- 
rect coupling to the thermal bath we shall consider the 
interaction of magnons with conduction electrons in a 
ferromagnetic metal. From a microscopic point of view 
the most probable process is the confluence of a magnon 
with wave vector k = and energy Uloq and a conduc- 
tion electron with wave vector ki 7^ and energy e(ki) 
into a conduction electron with wave vector k2 and en- 
ergy e(k2) = e(ki) -I- Tiloq. It is obvious that this process 
is forbidden for an ideal crystal, where the momentum 
conservation is valid \^2 = ki -I- 0. Such confluence pro- 
cesses exist only for non-uniform magnetic modes and 
demonstrate a strong fc-dependence for the relaxation 
rate (cx fc^) ||, |4|. 

The abovementioned confluence process can occur in 
the presence of defects, impurities or fluctuations which 
permit violation of momentum conservation in the crystal 
(Fig. la). This opportunity was qualitatively discussed by 
Kambersky and Patton Because of its importance, 
we shall consider this process in detail (see, Appendix 
B). The relaxation rate for the process can be written as: 



1V2 



duj, 



ek 



El/. 



(23) 

Wek)- 



This relaxation rate is linear with the defect concentra- 
tion Cdef and frequency dependent: it contains trans- 
formation terms (m^ -|- v'^)ujo = "l{Hxo + Hy^)/2. The 
temperature dependence can occur through the term: 
X^k (~^k/9wek)- An estimate of Eq.(p3|) for /kk' = 
/ = const and T <C ep/kB ^ 3 • 10"* K gives 



Y^7(^^a + ^yoJ ^4 ' 



(24) 



where m is the conduction electron mass and e^? is the 
Fermi energy. 



IV. RELAXATION VIA AN INTERMEDIATE 
DAMPED DYNAMIC SYSTEM 

There is a whole class of relaxation mechanisms that 
can not be analyzed as elementary processes. The energy 
loss in this case occurs via intermediate damped dynamic 
system. Here we shall consider the magnetization damp- 
ing via so-called, 'slow relaxing' impurities (see. Fig. lb 
and Refs. H], ||^, [|||, (29)). In this mechanism the 
magnetization motion modulates the impurity splitting 
(levels, see. Fig. lb). Thus the thermal equilibrium popu- 
lation of the energy levels varies and transitions between 
the levels occur (arrows in Fig. lb). There is a delay for 
these transitions due to a finite impurity relaxation time. 
This delay results in a magnetization oscillation energy 
loss. 



A general Hamiltonian describing the i mpu rity level 
modulation is given in Appendix C, Eq. ( |C2| ) . In the 
case of only the coherent motion of the host spins we 
have S(Rj -I- r^) = -(Fo/?i7)M, where Vq = V/N is the 
volume of elementary cell. Let us consider a two-level 
paramagnetic impurity with the energy 



(25) 



where floj is the splitting frequency, nj — Szj(Rj) -1-1/2 
is the upper lever po pula tion and j is the impurity index. 
Applying Eqs.(^) to (02), we can write the impurity level 
modulation as 



(26) 



where 



E 



Hy 



1/4 



Hx 
IT 



ya 



1/4 



(27) 



The kinetics of the impurity population is defined by 
the following equation: 



duj/dt — — — nT{i~lj)]. 



(28) 



Here Fy j is the impurity relaxation rate and nxi^j) = 
[eiip{hflj /ksT) + 1]^^ is the equilibrium population at 
frequency flj = figj + SQj(t). Taking into account 



nj{t) — riTi^oj) + 5nj{t), 



(29) 



we can solve Eq.(G8) and obtain 



<5n,(i)=F|| 



driT 



^jc(^) 
'11,7 - 



Fii 



(30) 



The dynamic equation for the normal mode interact- 
ing with impurities is obtained by substituting ( ^ ) into 
(pq). This yields an additional term to the basic mode 
Hamiltonian (^. Similar to (^, the dynamic equation is 
now: 



dc/dt ~ —iujQC — id{T-Limp,j /K) / dc* 
— —iajQC — i ^*Snj{t). 



(31) 



Substituting Eq.(|^) into (||) and (|3^), we obtain the 
equation of damped harmonic oscillator ||30[| : 

dc/dt + TjsrC = —i{uJo + ALL!sr)c + f{t). (32) 

Here the frequency shift is equal to 

F2, 



-El*. 



2 ( 9nr(Slo^j) 



(33) 



4 



and the relaxation rate is: 

,2 / 9nT(rioj)\ '^ory, 



da 



where 



1$, 



MsVo 



(34) 



Vo 



1/2 



To simpHfy (p^), we assume that rJo.j — ^^o a-nd Tn j = 
for aU impurities. Summing over j gives an average 
anisotropic exchange B: 



-XQ 

^ . (35) 



E 

i 

E 



y^^Zj,a:o(Rj,l'^ 

1/ 



(36) 



H 2 



where Nimp is the total number of impurities in the sam- 
ple and Zimp is the average number of magnetic neighbors 
for one impurity. Finally the relaxation rate ( ^^ is: 



^ AzimpB) exp(?if2o/A:Br) 

^sr — CimpO 



(37) 



^ifcsT [exp(ni7o/fcBr) + 1]^ 

7(g.o + g,jr|| 

^ p2 I , ,2 ' 

where c^mp = Ni^^p/N is the impurity concentration and 
S = MsVo/?i7 is the value of the host spin. This relax- 
ation rate exhibits a definite temperature dependance. 
As discussed in the previous section for magnon-electron 
scattering processes, the impurity relaxation rate also is 
frequency dependent, via the term j{Hxg + Hy^) (with 
identical dependence as in magnon-electron scattering). 



V. MACROSCOPIC DYNAMIC EQUATIONS 

A. Tensor form of magnetization damping 

We have seen that both direct and indirect magne- 
tization interactions with a thermal bath yield damped 
harmonic oscillator dynamics in the form: 



dc/dt + rjc — —i{LUo + Aw)c + f{t). 



(38) 



As discussed in g, ||], utilizing back transformation (|^), 
from (|3^) we can derive linearized equations for the trans- 
verse magnetization components: 



dt \ M 



yo 



-77 
iHxo 



-f] 
hyoit) 



Vo 



(39) 



Here 

hyoit) 
hxoit) 



\2-fMsV 
■ ( ^ ' 



tt" [nt) + f{t)], (40) 



1/4 



are the independent random fields. 

The nonlinear macroscopic dynamic equation for an 
anisotropic magnetic system which conserves the length 
of magnetization (|M| = M^) has the form ||], [0: 

^ = -7M X [Hoff + h{t)] - 7^ X [a -(M X Hcff)]. 



(41) 



Here a is a dimensionless damping tensor, which contains 
all necessary information about the system symmetry. In 
the vicinity of equilibrium one has 



^ 77 
a— — 

7 



l/Hy, 

l/H^o 




(42) 



and Eq.(^) is equivalent to (g9|). 

Note that the tensor form of damping in (^Tj) is a nat- 
ural description of an anisotropic system. It reflects an 
elliptical motion of magnetization and contains just one 
relaxation parameter 77. The a should not be confused 
with a damping tensor which contains two or more inde- 
pendent relaxation parameters. 



B. Linearized Landau-Lifshitz equation 

The Landau-Lifshitz equation with random fields h(t) 
has the form: 

— = -7M X [Heff + h{t)] - j^M X (M X Heff) , 

(43) 

where a is a dimensionless damping parameter. Us- 
ing (|^), we can calculate the effective field Hcff = 
—d{£/V)/d'M. and write down the linearized equations 
for the transverse magnetization components {Mz^ ~ 
M,): 



d { 



dt \ M, 



ya 



-a-fHx„ -jHy„ \ / j ^^^^ 
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We see that the non-diagonal terms in (^J) and (39), 
as expected, coincide with each other, respectively. The 
diagonal terms, responsible for relaxation, in general, are 
different (-ff^o 7^ Hy^). These damping terms are equal 
only in the special case when Hx^ — Hy^ (see. Appendix 
A). 

The Eq. (H) can be rewritten in terms of normal modes 
using (0) as 



where from (@) 



dc/dt + ric + rjc* = —iloqc + /(t). 



(45) 



where 



il = a^{Hx,+Hy,)/2, ?j=a^{Hx,-Hy„)/2. (46) 

We see that Eq. ( p5| ) contains one extra fast oscillating 
term 77c*. As shown in Ref. |P, this term changes the 
power spectral shape. However, in the case when the 
thermal term is negligible, the fast oscillation term 77c* 
effectively vanishes by averaging. In those cases the LLG 
dynamics can be used as an approximation. 

The most important influence of this parasitic fast os- 
cillation appears in the case when we study the effect 
of a random force f{t) in (p5|); this term gives effec- 
tively higher magnetization noise in the system. ^ , . 
The thermal bath produces both magnetization relax- 
ation and thermal fluctuations whose level is defined by 
fluctuation-dissipation theorem. Both damping and ther- 
mal noise appear as a result of microscopic spin fluctu- 
ations, therefore a connection of phenomenological dy- 
namic magnetization equations with microscopic physics 
is necessary. We can quote van Kampen |3l|| : "The moral 
is that one cannot give a stochastic description of a physi- 
cal system simply by taking one of the familiar equations 
for the sure case and declaring some of the coefficients 
to be random. One has to start from the fundamental 
equations..." . 



VI. ANALYSIS OF FMR EXPERIMENTAL DATA 

FMR linewidth measurements have been widely uti- 
lized to determine the relaxation rate for low-level linear 
excitations |32| . In this section we will use the relaxation 
rates for microscopic mechanisms analyzed in Sees. Ill 
and IV to analyze FMR data in soft thin films (NiFe, 
permalloy). 

For small relaxation compared to the resonant fre- 
quency the linewidth is given by (e.g., 0): 



(47) 



for all microscopic processes. Usually in FMR experi- 
ments, the field swept linewidth Ai? is measured. There 
is a simple relation between A_ff and frequency linewidth 
Aw: 



AiJ = 



Aw 



duo/dHo ' 



(48) 



duo/dHo - -f\Hx„ + Hy,)/2ujQ, 



(49) 



which is valid for any direction of the applied magnetic 
field, including those parallel and perpendicular to the 
film plane. 

Let us summarize the linewidths obtained in Sections 
III and IV. Using Eqs. (||), (|^, (g) and (liS 
obtain the field swept linewidths: 

magnon- electron confluence process 



Cdef ^0 m^Vplfl ep 
47r3 T fi-* 



(50) 



slow relaxing impurity process 



AH, 



hj ksT [e^p{hno/kBT) + 1]2 F^ 



(51) 



Note that the linear frequency dependence for AHm-e 
was predicted in Ref. ||25[| . 

The temperature dependence of AH in permalloy thin 
films with the static magnetic field in the film plane has 
been measured by Patton and Wilts As seen in 

Fig. 2, AH exhibits a strong frequency dependence with 
a maximum in the vicinity of Tmax ~ 80 K. The tem- 
perature maximum shifts to slightly higher temperatures 
with increasing FMR frequency. Such a non-monotonic 
temperature dependence of AH is typical for slow relax- 
ing impurities (see, @). However, the slow- 
relaxing impurity process alone can not describe the ex- 
periment . One can assume that the magnon-electron 
confluence process plays a role of a relaxation rate "base- 
line" . Thus for analysis we combine the slow-relaxing 
impurity and magnon-electron confluence processes to- 
gether: 

AH = AHsr{T,UJo,Cimp) + AHm-eil^O,Cdef)- (52) 

According to Eq. (jsj) , the temperature T^ax of linewidth 
maximum is given by 



ril (Tmax) — 1^0 • 



(53) 



The impurity damping Fn is expected to increase mono- 
tonically with temperature [ p9[ . Thus Tmax will increase 
with increasing luq in agreement with experiment. The 
solid lines in Fig. 2 represent a theoretical fit with Eq. (p2[) . 
We see that the temperature (and frequency) dependence 
( p^ ) describe well the range from about the peak to about 
room temperature. 

For the slow relaxation mechanism the fit gives 
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AiJ,^[Oe] 



1.9- 10" 



exp 



( 100. ^ 

\Tm) 



"0 

Til 



r[K] 



(54) 



where typical parameters for slow-relaxing impurity were 
used: Sf^o/fcs = 100 K, Zi„,pB/kB = 25.3 K, 5* ~ 1 and 
Cimp — 10~^. The spin-lattice relaxation rate variation 
with temperature was assumed to be r||(T) = cyT^ with 
a fit of C||/27r = 0.8-10"^ GHz/K^. The impurity interac- 
tions with conduction electrons are likely to be responsi- 
ble for such a temperature dependence. The spin-lattice 
relaxation time at room temperature l/r||(300 K) ~ 2.2 
psec agrees with typical experimental data p9| . 

For the magnon-electron confluence process the fit 
gives 



AHm-e[Oc\ ~ 2.1 (cjo/27r[GHz]) , 



(55) 



Substituting m ~ 10 ^ g, Vo — 10 cm'^ and 
7 — 1.76 10^ s~^Oe"^ into Eq.(pO[) and comparing with 
Eq.dl), we obtain Cdef\fif?^F ^ 7 • lO^^s erg^. Tak- 
ing defect concentration Cdef — 10^'^ and Fermi energy 
= 3 eV, one gets an estimate for the interaction am- 
plitude hf ~ 2.4 eV, which is typical for ferromagnetic 
metals ||34| ]. 

Patton et al [|l2| have studied the frequency depen- 
dence of AH in thin permalloy films (17-48 nm) at room 
temperature (see, Fig. 3). They have found that the field 
linewidth has a linear frequency dependence AH cx loq 
for the case when the external magnetic field is parallel 
to the film plane (open circles in Fig. 2). For the perpen- 
dicular case (black circles) the linear dependence is valid 
just for high FMR frequencies uo/^n > 8 GHz and it is 
saturated at smaller frequencies. Recent data show 
similar results except that for the field perpendicular to 
the film plane the saturation region occurs only at ex- 
tremely small frequencies (Li;o/27r < 2 GHz). Essentially 
for a wide range of frequencies the field swept linewidths 
are identical for the two field orientations and strictly lin- 
ear with resonant frequency. It is clear that this linearity 
can be explained by the magnon-electron confluence pro- 
cess (|55|). At room temperature the impurity relaxation 
rate is r||/27r ~ 72 GHz so that this mechanism yields a 
linewidth (^4|) that is also linear with frequency at least 
to about 50 GHz. Taking (||), (||) and (||) at T = 300 
K we obtain: 



AH[Oe] ~ 2.3 (wo/27r[GHz]) . 



(56) 



This linear dependence is shown in Fig. 3 by a solid line 
and demonstrates good agreement with experiment. 

It is interesting that Eq.(^2|) can be also used for in- 
terpretation of experimental field linewidth in ultra-thin 
permalloy films | |36[ | . The thickness dependence at room 
temperature for two different FMR frequencies is shown 
in Fig. 4. Assuming that the surface magnetic atoms be- 
haves as impurities, the concentration Cimp increases with 



decreasing film thickness as Cimp oc do/d, where do ~ 3 A 
is the linear size of the single-atomic magnetic layer. The 
theoretical thickness dependence for typical microscopic 
parameters (as above) can be written as 

AiJ[Oe] = (1.3 + 53/d[A]) (cjo/27r[GHz]) . (57) 



The solid lines in Fig. 4 by plot Eq.(|57|) and demonstrate 
good agreement with experimental data. This is an alter- 
native interpretation of [^6| , where two-magnon scatter- 
ing on surface defects (with thickness dependence l/d^) 
|2l| plus a baseline thickness independent process were 
suggested. In order to evaluate these two interpretations, 
FMR linewidth should be measured at lower tempera- 
tures, where the role of slow-relaxing impurities increases. 



A. Comment: data fitting by the LLG equation 

The Landau-Lifshitz-Gilbert equations were originally 
introduced (a) for small magnetization motions (small 
deviations from the equilibrium) and (b) for the case 
of high magnetic symmetry (axial symmetry) with an 
isotropic damping fitting parameter a (or "damping con- 
stant" ) ||l[ , . Nevertheless both "a" and "b" conditions 
are usually violated; these equations have been used for 
large magnetization motions such as magnetization re- 
versal and for a description of magnetization dynamics 
in systems with lower magnetic symmetry (e.g., in thin 
films) . 

In the LLG phenomenology the linewidth is given by 
(e.g., §): 



Auj = a-/ {H,,+Hy,). 



(58) 



From Eqs. (^) and ( p8| ) we find that the corresponding 
field swept linewidth is: 



AH = aujo/^. 



(59) 



This formula has been used to obtain the phenomenolog- 
ical damping parameter a from experimental data AH 
and frequency loq. 

Assuming a constant a leads, according to (p9[), to a 
linear frequency dependence for AH, in agreement with 
the room temperature data shown in Fig. 3. This agree- 
ment has been used to claim the validity of the LLG equa- 
tions even for an anisotropic system, e.g. a thin film [ p5[ . 
This dependence is obviously coincidental. The LLG ap- 
proach is purely phenomenological and can be related to 
physical processes only in the case of axial symmetry. As 
we have shown in this paper, detailed microscopic analy- 
sis yields dynamic equations which depend on the system 
symmetry. Microscopic analysis with appropriate relax- 
ation mechanisms explain a wide variety of experimental 
data and, as well, can be used to predict dynamic mag- 
netization phenomena. 
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VII. DISCUSSION 

In this paper we have examined the effect of general 
forms of microscopic relaxation mechanisms on magneti- 
zation dynamics. Only the case of linear relaxation has 
been considered when the deviations of magnetization 
from equilibrium are relatively small. We have shown 
for two broad classes of relaxation mechanisms that dy- 
namic relaxation may be expressed in terms of a damped 
harmonic oscillator equation driven by thermal fluctua- 
tions. The coordinates of the harmonic oscillator reflect 
a transformation of the transverse magnetization coor- 
dinates to circular motion. The transformation directly 
reflects the magnetic anisotropy. The relaxation time or 
damping parameter contains details of the microscopic 
relaxation mechanisms and also, in general, information, 
via the transformation, about the magnetic anisotropy. It 
is important to note that by anisotropy we do not mean 
the for of the interaction with the thermal bath. Here 
the anisotropy that yields a tensor form is macroscopic 
and arises, for example, from the sample shape (e.g., a 
thin film) or the crystalline symmetry. 

Specific mechanisms of slow relaxing impurities and 
magnon-electron coupling were analyzed. It was shown 
that a sum of these two effects explains very well FMR 
field swept linewidth measurements versus temperature, 
frequency and film thickness. The fitting parameters 
agree well with previous estimates. 

Transformation of the harmonic oscillator back to mag- 
netization components results in a macroscopic form that 
is different from that of LLC The damping term, al- 
though proportional to the relaxation time or damping 
parameter, depends on the magnetization and effective 
field through a tensor, whose components reflect the mag- 
netic anisotropy. LLG results only for the case of uniaxial 
system anisotropy about equilibrium. The LLG equa- 
tions for a fixed damping parameter a do show a linear 
frequency dependence of the linewidth, in agreement with 
data. However, this is just an artifact: the LLG damping 
can not be derived in terms of a microscopic relaxation 
rate. 

In 10] and earlier in it was shown that a magnetic 
system capable of non-uniform excitations about (a pos- 
sibly) non-uniform equilibrium also resulted in damped 
harmonic oscillators. In this general case the eigen- 
modes of the general gyromagnetic motion tensor (for 
a discretized system) yield independent harmonic oscil- 
lators with corresponding resonance frequencies, relax- 
ation rates and stochastic thermal energy. Examples 
were given for a thin film sensor in a GMR transducer 
and in [ p8| experimental fitting showed a better agree- 
ment than with simply adding energy to each discretized 
computation cell and applying LLG with random fields 
to each cell. Note that the dynamic normal mode ap- 
proach with phenomenologically introduced damping (as 
a small imaginary addition to the frequency) has been 
successfully used in the theory of parametrically excited 



spin waves (see, e.g., ||, [|o|, g). 

We have examined microscopic magnetization dynam- 
ics for two broad classes of interactions with a thermal 
bath. Further work will study the role of anisotropy in 
the case of magnetic multilayers and other indirect 
processes M. 
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APPENDIX A: ELLIPTIC MAGNETIZATION 
MOTION 



Consider the following energy density: 

£p/V = K\\ sin^ e + K^{Mx/Msf - M • Hq. (Al) 

Here K\\ > and K± > are the "easy" {z) and "hard" 
(x) uniaxial anisotropy constants, 9 is the angle between 
the magnetization M ={Mx, My, Mz) and the easy axis 
of anisotropy. Hq = Ho{0, sinOn, cos9h) is the exter- 
nal magnetic field. The equilibrium magnetization state 
6* = 6*0 is defined by the condition d£p/d9 = 0, or. 



i?^" sin 26*0 = 2Hosm{0H - 0q), 



(A2) 



where H^^'^ = 2K\\/Ms. In order to describe small mag- 
netization oscillations it is necessary to transform the co- 
ordinate system and magnetization as follows 

y = 2/0 cos 6*0 + zq sin 6*0, (A3) 
z — —yo sin 6*0 + zq cos Bq, x = xq, 

My = My„ cos 00 + Mzg slu 6*0 , ( A4) 

Mz = -My, sin 00 + Mz„ cos 0o , = . 



and 



The equilibrium state in these coor- 
dinates is M=(0, 0, Ms). Taking into account trans- 
formation (|a3) and relations = - - 
and Mz, ~ Ms - (Mj^ + M2J/2M,, the energy (^) can 
be written in the form (ih, where 



Hx, = i?^" cos2 00 + hI^^ + Ha cos{0h - 0o), 



Hy„ ^ hI^^ COS20O - 



Hocos{0H - 0o)- 



(A5) 



and H£' = 2Kjl/Ms. From (|AJ) and (|AJ) it follows 
that = Hyg only in two special cases: 1) the case 
of spherical symmetry and 2) the case of uniaxial sym- 
metry when the external magnetic field and equilibrium 
magnetization are oriented along the easy axis. 



8 



APPENDIX B: MAGNON-ELECTRON 
SCATTERING WITHOUT MOMENTUM 
CONSERVATION 

The interaction of the uniform magnetic precession (a, 
a*) with electrons is assumed to be of the form: 



where 



Hi = {D^a + Da*) /Vn, 



(Bl) 



(B2) 



k,k' j 



Here /kk'(rj) = |/kk'l exp[i(/)(rj)] describes the scattering 
process in the vicinity of crystal defect at the point Vj, 
(j>{rj) is the phase, dj^. and c?k are the fermion operators 
describing creation and annihilation of electron with wave 
number k and frequency LUek- 

Using transformation (^) to a normal mode, we can 
rewrite Tii as: 



k,k' J 



APPENDIX C: INTERACTIONS WITH 
IMPURITIES 

Let us consider two-level impurities as effective spins 
Sj = 1/2. The anisotropic exchange Hamiltonian be- 
tween the impurities and the neighboring host spins can 
be written in a general form: 

^ex - ^^Ba,,a(Rj,r.)Sa,(Rj)^a(Rj +r,). (CI) 

j.v aj^a 

Here Baj,a(Rj , Ti.) 

are the anisotropic exchange integrals, aj = Xj , yj , Zj 
are the local principal axes for the j-th impurity "spin" 
Sj(Rj) = {sxj (Rj), Sy^ (Rj), (Rj )} located at R^-, and 
S (R, + ) = { (R, + ) , Sy„ (R, + I"!/), Szq (R,+r,)} 
are the components of the host spin located in the vicin- 
ity of j-th impurity. 

The Hamiltonian describing the impurity level modu- 
lation follows from Eq.(Cl): 



'^ex,slow — T ^ ^ \Bzj ,xo(^j ,'^v)Sxo(^i + 1",^) (C2) 



(B3) 



+-B;^j,ao {^o,T^y)Syo (Rj + r^)] s^^ (R^). 



where 



*kk'(rj) = it/kk'(i"j) + t'/k'k(i"j) 



(B4) 



is the amplitude of magnon-electron scattering. 

The relaxation rate is defined by the golden Fermi rule 

as 



P X! E l*k'kP("-k - "k')'5(^ek' - t^ek - ^o) 



N3 



k,k' 



^^^aEEKiz-i^+^'i/-'!^ 

k,k' j 

+™|/k'k/kk'| cos[20(rj)]} 
x(nk - ni^')SiuJe-k' - t^ek - t^o)- 



(B5) 



Here TTk = (c^k'^k) is the Fermi occupation number If we 
assume that the scattering phases (/>(rj) are random and 
|/k'k| = l/kk'li then the uv term vanishes and Eq.(B5) 
becomes 



7V2 



El/. 



kk' I 



fc,fc' 



(B6) 



Here Cdef = Ndef/N, Ndef is the total number of de- 
fects. After the following simplifications: — rik' — 
—ujodfik/dojeu and ujq <§; luq, we obtain Eq. (p3|) . 
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Figure captions 

Fig.l Relaxation processes: a) confluence of magnon 
and conduction electron into conduction electron in the 
presence of defect, b) 'slovif-relaxation' mechanism. via 
modulation of impurity levels by magnetization oscilla- 
tion. 

Fig. 2 The linewidth temperature dependence for three 
FMR frequencies. Points are experimental data of Ref. 
|ll|. Solid lines demonstrate theoretical fit with (54) and 



Fig. 3 FMR frequency dependence at room tempera- 
ture. Open circles correspond to parallel and black point 
to perpendicular field orientations . Solid line demon- 
strate theoretical fit (|5^). 

Fig. 4 The linewidth thickness dependence for two FMR 
frequencies at room t emp erature . Solid lines demon- 
strate theoretical fit (|57|). 
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